AD-A041  166  NAVAL  POSTGRADUATE;  SCHOOL  MONTEREY  CALIF 

THE  USE  OF  DMSP  AND  SMS-2  DIGITAL  SATELLITE  DATA  FOR 
MAR  77  C B IHLI>  R J RENARD 
NPS-63IH77031 


F/G  4/2 

IDENTIFYIN—ETC(U) 

NL 


UNCLASSIFIED 


DDE  FILE  COPY,  AO  A 0 4 1 1 6 6 


NPS-63Ih77031 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 


THE  USE  OF  DMSP  AND  SMS-2  DIGITAL 
SATELLITE  DATA  FOR  IDENTIFYING  MARINE 
FOG  IN  THE 

EASTERN  NORTH  PACIFIC  OCEAN  AREA 


Carl  Bolton  Ihli,  Jr, 


March  1977 


Thesis  Advisor: 


R . J . Renard 


Approved  for  public  release;  distribution  unlimited. 


SECURITY  CLASSIFICATION  of  This  FACE  (Whm%  Oau 


REFORT  DOCUMENTATION  PAGE 


NPS-63Xh77031 


. TITLE  (And  SuOtltlo)  ' 

The  Use  of  DMSP  and  SMS-2  Digital  Satel- 
lite Data  for  Identifying  Marine  Fog  in 
the  Eastern  North  Pacific  Ocean  Area 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  RECIPIENT'S  CATALOG  MUMICR 


■ TYPE  OF  REPORT  A PERIOD  COVERED 

Master's  Thesis; 

March  1977 


. PERFORMING  ORC.  REPORT  NUMBER 


7.  authors; 

Carl  B.  Ihli,  Jr. 

In  conjunction  with 
Robert  J.  Renard 


*.  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Naval  Postgraduate  School 
Monterey,  California  93940 


II.  CONTROLLING  OFFICE  NAME  ANO  ADORES* 

Naval  Air  Systems  Command,  370C 
Washington,  D.  C.  20360 


12.  REPORT  DATE 

March  1977 


IS.  numeer  of  pages 

99 


1*.  MONITORING  AGENCY  NAME  a AOORESSflf  Ollloronl  from  Controlling  Ollleo)  | IS.  SECURITY  CLASS.  (»l  Mila  rOron) 


Naval  Postgraduate  School 
Monterey,  California  93940 


Unclassified 


<Sa.  DECLASSIFICATION/  OOPNGR  AOlNG 
SCHEDULE 


I*.  DISTRIBUTION  STATEMENT  (ol  UUa  Raparl J 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  Mia  aAalraal  omtmroO  M Hook  70,  II  dllfarail  hom  Moron) 


t*.  KEY  WOROS  (Contitnto  on  rororoo  mlOo  II  nocoooorr  and  loom  Mr  Or  Hook  mA«l 

f og ; marine  fog;  weather  satellite;  North  Pacific  Ocean;  ocean 
fog;  DMSP  satellite  data;  SMS-2  satellite  data;  geosynchronous 
satellite;  meteorological  satellite;  marine  visibility. 


20.  ABSTRACT  (Conilmuo  an  rororoo  alda  II  notoooorr  and  I dam  Iff  Or  AlaaA 

’Defense  Meteorological  Satellite  Program  (DMSP)  and 
Synchronous  Meteorological  Satellite  (SMS-2)  digital  satel- 
lite data,  both  visual  and  infrared,  are  processed  and 
analyzed  in  context  with  transient/stationary  synoptic  ship 
reports,  in  an  effort  to  identify  areas  of  marine  fog  occur- 
rence. The  eastern  North  Pacific  Ocean  was  chosen  as  the  — 


hex 


50  i j a"  7i  1473 


COITION  Of  1 NOV  ••  It  OOtOLKTK 
S/N  0102-014-  6401  I 


tecum  TV  Cl  asaipicatiom  or  this  m«o  irm*,  n.,.  CmeMi 


cr^X 


■area  of  study;  for  v 
to  5-9  August  1976. 
of  visual  and  infrar 
of  five  categories  o 
ranges  of  visual  and 
ciated  with  each  of 
in  the  study  area  ar 
of  the  cloud  base  is 
data  points,  based  o 
physics  and  empirici 
Wisconsin,  Insofar 
identifying  marine  f 
vations  is  concerned 
sidered  suggestive  b 
limited  data  sample, 
work  at  the  Naval  Po 
orbiting  satellite  d 


arious  reasons  the  period  was  limited 
Histograms  and  two-dimensional  plots 
ed  data  are  generated  as  a function 
f fog/no-fog  occurrence.  Optimal 
infrared  digital  count  values  asso- 
the  five  categories  of  fog  occurrence 
e proposed.  In  addition,  the  height 
computed  for  all  processed  SMS-2 
n a formulation  derived  from  radiation 


sm 

available 

fr 

om 

the 

Universi 

ty 

of 

as 

advanc i 

ng 

th 

e 

skill 

of  obje 

c ti 

ve  ly 

og 

areas  f 

r om 

wea 

ther 

satellit 

e obser 

.,  the  resu 

Its 

o 

f 

this 

study  ar 

e c 

on- 

IUt 

not  def 

ini 

ti 

ve 

, in 

part  due 

to 

the 

The  rela 

tio 

n 

of 

the 

study  to 

pr 

ior 

'Stg 

raduate 

Sc 

ho 

ol 

, us  i 

ng  NOAA 

pol 

ar- 

,ata 

, i s no 

ted 

• 

/V 


DD  Form  14i3 
1 Jan  73 

S/N  0102-014-6601 


*CCU«lTV  CLAMIPICATiOn  OP  TMI*  PiMWw  OAfP 


r 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


Rear  Admiral  Isham  Linder  Jack  R.  Borsting 

Superintendent  Provost 


This  thesis  prepared  in  conjunction  with  research 
supported  by  the  Naval  Air  Systems  Command  370C. 

Reproduction  of  all  or  part  of  this  report  is 
authorized. 


Released  as  a 
Technical  Report  by: 


R.  R.  Fossum 
Dean  of  Research 


Approved  for  public  release;  distribution  unlimited 


The  Use  of  DMSP  and  SMS-2  Digital 
Satellite  Data  for  Identifying  Marine 
Fog  in  the 

Eastern  North  Pacific  Ocean  Area 


by 


Carl  Bolton  Ihli,  Jr. 
Lieutenant,  United  States  Navy 
B.S.,  United  States  Naval  Academy,  1968 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  METEOROLOGY  AND  OCEANOGRAPHY 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
March  1977 


Author 


Approved  by: 


JCtESMOU  t«  ! 

•IIS 

»«'!*  Seel ili 

Y 

aoc 

n 

KltMSOIIHCE* 

□ 

jus  i if  it;  no* 

'3IHTT  60  0 0 
,r  rPECIH 


J //AA4t*st/  ' 

Thesis  Advisor 


3 


n n - 


f7V- 


\u< 


hi,  r 

UUiiic 


ABSTRACT 


Defense  Meteorological  Satellite  Program  (DMSP)  and 
Synchronous  Meteorological  Satellite  (SMS-2)  digital  satel- 
lite data,  both  visual  and  infrared,  are  processed  and 
analyzed  in  context  with  transient/stationary  synoptic  ship 
reports,  in  an  effort  to  identify  areas  of  marine  fog  occur- 
rence. The  eastern  North  Pacific  Ocean  was  chosen  as  the 
area  of  study;  for  various  reasons  the  period  was  limited 
to  5-9  August  1976.  Histograms  and  two-dimensional  plots 
of  visual  and  infrared  data  are  generated  as  a function  of 
five  categories  of  fog/no-fog  occurrence.  Optimal  ranges 
of  visual  and  infrared  digital  count  values  associated  with 
each  of  the  five  categories  of  fog  occurrence  in  the  study 
area  are  proposed.  In  addition,  the  height  of  the  cloud 
base  is  computed  for  all  processed  SMS-2  data  points,  based 
on  a formulation  derived  from  radiation  physics  and  empiri- 
cism available  from  the  University  of  Wisconsin.  Insofar 
as  advancing  the  skill  of  objectively  identifying  marine 
fog  areas  from  weather  satellite  observations  is  concerned, 
the  results  of  this  study  are  considered  suggestive  but  not 
definitive,  in  part  due  to  the  limited  data  sample.  The 
relation  of  the  study  to  prior  work  at  the  Naval  Postgraduate 
School,  using  NOAA  polar-orbiting  satellite  data,  is  noted. 
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I 


INTRODUCTION  AND  BACKGROUND 


The  occurrence  of  fog  is  a threat  to  all  ocean-going 
vessels,  civilian  and  military  alike.  In  naval  applica- 
tions, aircraft  operations  as  well  as  surface  ship  opera- 
tions are  affected  by  the  visibility.  A study  concerning 
materiel  dollar  cost  and,  more  importantly,  losses  of  human 
life  over  a five-year  period  further  emphasize  the  impact 
that  marine  fog  has  on  naval  operations  (Wheeler,  1974). 

Accurate  depiction  of  fog  regions  over  the  open  ocean 
on  a real-time  basis  would  be  of  significant  value  in  any 
naval  operation  in  which  horizontal  visibility  is  a limit- 
ing parameter.  Ship  routing,  by  way  of  example,  would  cer- 
tainly take  into  account  any  reliable  methods  of  fog  diag- 
nosis and  prediction  should  such  products  exist. 

Currently,  only  Fleet  Numerical  Weather  Central, 
Monterey,  California  (FNWC)  generates  objective  fog  fore- 
casts. Their  product  (FTER)  forecasts  fog  occurrences  for 
all  areas  of  the  Northern  Hemisphere.  FNWC  promulgates  the 
computerized  probability  of  fog  twice  daily  (0000  and  1200 
GMT)  for  forecast  intervals  up  to  72  hours.  These  forecasts 
are  based  on  a statistical  processing  of  certain  model  out- 
put parameters  related  to  fog  occurrence  (U.  S.  Naval 
Weather  Service  Command,  1975)  . 

The  use  of  the  meteorological  satellite  may  lead  to  an 
improvement  in  diagnosing  marine  fog  and  concomitant 
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increases  in  predictive  accuracy.  Research  conducted  at 
the  Naval  Postgraduate  School,  Monterey,  California  (NPS) 
since  1974,  made  use  of  visual  (VIS)  and  infrared  (IR) 
weather  satellite  imagery  from  polar  orbiting  satellites 
to  discern  the  presence  of  marine  fog.  Initial  results 
indicated  the  need  for  further  investigation  into  the  use 
of  digital  satellite  data  in  detecting  marine  fog  (Wallace 
and  Renard,  1975).  Continuing  investigations,  using  NOAA-2 
digitized  satellite  data,  established  optimal  ranges  of 
count  values  for  fog  occurrence  for  both  VIS  and  IR  digital 
data  (Hale  and  Renard,  1975).  However,  the  level  of  skill 
in  applying  these  findings  indicated  the  need  for  follow-on 
research. 
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II 


OBJECTIVES  AND  APPROACH 


From  the  research  cited  in  Section  I there  appeared  to 
be  two  directions  that  future  studies  could  take.  One 
possible  approach  involved  the  analysis  of  more  data  observed 
by  the  same  or  another  polar  orbiting  satellite.  To  this 
end  the  Defense  Meteorological  Satellite  Program  (DMSP) 
satellite  was  selected,  largely  due  to  its  application  for 
military  oriented  operations  in  general,  and  because  its 
output  will  be  eventually  readily  available  to  FNWC  for 
operational  purposes.  Although  the  DMSP  satellite  is  con- 
trolled by  the  United  States  Air  Force,  its  products  are 
in  wide  usage  by  Naval  Oceanography  and  Meteorology  Units 
as  evidenced  by  the  DMSP  satellite  installations  at  several 
land  and  afloat  sites. 

The  second  possible  approach  was  to  obtain  data  from 
a geosynchronous  weather  satellite  and  thereby  have  the 
additional  advantage  of  nearly  continuous  coverage  of  an 
ocean  area.  This  second  option  was  also  selected  and  the 
Pacific  Ocean  Synchronous  Meteorological  Satellite  located 
over  the  equatorial  Pacific  Ocean  area  (SMS-2)  was  chosen. 

Thus,  the  primary  objective  of  the  subject  research  is 
to  ascertain  the  value  and  anplicability  of  both  DMSP  and 
SMS-2  satellite  data  in  identifying  areas  of  marine  fog 
over  the  open  ocean. 
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ll 

The  approach  taken  in  pursuit  of  this  objective  was  to 
statistically  process  and  interpret  DMSP  and  SMS-2  digital 
satellite  data,  with  synoptic  weather  reports  from  transient 
and  stationary  ships  used  as  "ground  truth"  verification. 

The  immediate  goal  was  to  establish  a limited  range  of  VIS 
and/or  IR  digital  count  values  which  would  allow  fog/no-fog 
discrimination.  In  addition,  usable  statistical  correla- 
tions of  satellite  and  convent iona 1 ly- reported  weather  para- 
meters were  sought. 
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DATA 


A.  AREA  OF  STUDY 

The  previously  cited  marine  fog  research  conducted  at 
NPS  was  confined  to  the  North  Pacific  Ocean  area.  The  area 
selected  for  this  study,  as  dictated  by  the  SMS-2  satellite 
coverage,  was  the  eastern  portion  of  that  ocean,  with  the 
following  boundaries:  30-60N;  115W-180.  Although  this  area 

includes  a portion  of  the  west  coast  of  the  United  States, 
only  ocean  areas  were  considered  in  the  study  (see  Fig.  1). 

B.  TIME  PERIOD  FOR  STUDY 

The  limiting  factor  in  the  determination  of  a time 
period  of  investigation  proved  to  be  the  period  for  which 
the  digital  SMS-2  satellite  observations  could  be  obtained 
from  the  receiving  site  at  the  University  of  Wisconsin, 
Madison,  Wisconsin.  Initially,  based  on  climatology  and 
previous  research,  a period  during  the  month  of  July  was 
requested.  July,  climatologically , provides  the  highest 
frequency  of  fog  occurrence  over  the  North  Pacific  Ocean 
(Willms,  1975).  Unfortunately,  the  month  of  July  was  pre- 
empted by  another  user.  However,  the  period  1-15  August 
1976  was  available  and  arrangements  were  made  for  its  ac- 
quisition. However,  due  to  weather  problems  at  Madison, 
Wisconsin,  the  receiving  site  was  shifted  to  the  Cambridge 
Research  Laboratory,  Cambridge,  Massachusetts.  Some 


equipment  difficulties  at  the  east  coast  receiving  station 
eliminated  any  data  after  9 August  1976,  so  that  SMS-2  data 
were  available  for  the  study  only  in  the  1-9  August  1976 
time  frame. 

DMSP  data  presented  different  problems.  As  had  been 
advertised,  DMSP  data  were  made  available  from  the  United 
States  Air  Force  Global  Weather  Central,  Offutt  Air  Force 
Base,  Nebraska,  during  the  proposed  1-15  August  time  period. 
Unfortunately,  the  data  as  received  at  NPS  did  not  cover  the 
entire  15-day  period.  On  some  days  data  were  unavailable; 
on  other  days  either  IR  or  VIS  data,  but  not  both,  were 
available.  On  certain  days  when  both  VIS  and  IR  data  were 
available,  the  data  in  the  study  area  were,  in  part,  from 
various  observation  times  differing  by  12  or  more  hours  or 
from  the  NOAA-4  satellite  if  DMSP  data  were  missing. 

Considering  the  limitations  imposed  on  the  study  by  both 
satellite  systems  and,  after  examining  all  data  once  re- 
ceived at  NPS,  a study  period  of  only  five  days,  5-9  August 
1976,  was  possible. 

C.  SYNOPTIC  WEATHER  REPORTS  FROM 

TRANS I ENT/ S TATI ONA  RY  SHIPS 

Primary-time  synoptic  weather  observations  for  5-9 
August  1976  were  provided  by  FNWC  in  the  main,  with  several 
missing  synoptic  times  data  available  through  archived  data 
sources  at  the  National  Climatic  Center  (NCC) , via  the 
co-located  Naval  Weather  Service  Detachment  (NWSD) , at 
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Asheville,  North  Carolina.  Surface  ship  synoptic  reports 
established  the  basic  "ground- truth"  for  the  existence  or 
non-existence  of  marine  fog. 

Initially,  the  four  primary  synoptic  times  of  0000, 

0600,  1200  and  1800  GMT  were  considered.  Later  in  the 

study  it  became  necessary,  due  to  other  limitations,  to 
consider  only  0000  and  1800  GMT  reports.  The  reason  for 
this  further  limitation  will  be  made  apparent  in  a subse- 
quent section. 

Of  the  information  available  in  a standard  surface  ship 
report  all  data  were  initially  recorded  in  tabular  format 
(see  Fig.  2,  for  example).  After  filtering  out  reports 
that  were  outside  the  geographic  limits  of  consideration, 
a total  of  1335  ship  reports  were  available  at  0000,  0600, 
1200  and  1800  GMT  for  the  five-day  study  period.  After  the 
synoptic  times  0600  and  1200  GMT  were  removed,  a total  of 
825  ship  reports  were  available  for  analysis. 

D.  DMSP  SATELLITE  DATA 

DMSP  digital  data  on  magnetic  tape  were  furnished  by 
the  Air  Force  Global  Weather  Central  (AFGWC)  at  Offutt  Air 
Force  Base,  Nebraska,  through  the  auspices  of  the  Naval 
Weather  Service  Detachment  located  there.  Data  as  stored 
on  the  magnetic  tape  were  depicted  on  the  standard  I-J 
grid,  as  used  by  FNWC.  This  grid  is  superimposed  on  a polar 
stereographic  projection  of  the  earth,  true  at  60  N latitude, 
with  the  pole  at  its  center.  It  consists  of  63  columns 
, and  63  rows  of  equally  spaced  points.  The  I-J  grid  coverage 
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of  the  study  area  is  included  in  Fig.  3.  The  distance  be- 
tween points  on  the  grid  is  related  to  distance  on  the  map 
by  the  map  scale  factor:  1 + sin  60°/l  + sin  (latitude). 

This  distance  at  60  N is  206  nmi  and  reduces  to  110  nmi  at 
the  equator. 

There  are  a total  of  64  by  64  DMSP  data  points  (i.e. 
digital  count  values)  within  one  single  square  of  the  I-J 
grid  (Fig.  3)  for  a total  of  4096  data  points.  The  study 
area  included  153  whole  or  partial  I-J  grid  squares  so  that 
the  total  number  of  data  points  available  at  any  one  report- 
ing time  over  the  study  area  was  very  close  to  600,000. 

Fig.  4 is  a sample  64  x 64  grid  of  two-digit  digital  count 
values  (DCV's) . Table  I summarizes  the  spacing  of  data 
points  (data  interval)  within  a unit  I-J  grid  square  as  a 
function  of  latitude. 

The  digital  data  that  appears  at  each  of  the  data  points 
in  the  64  by  64  matrix  ranges  from  DCV's  of  0-63  indicative 
of  64  shades  of  grey  for  both  visual  (VIS)  and  infrared  (IR) 
digital  data.  For  VIS  data,  the  lower  DCV's  represent  areas 
with  a very  low  albedo,  such  as  the  sea  surface.  In  typical 
DMSP  VIS  imagery  these  surfaces  appear  very  dark  grey  or 
black.  Dense  cloud  layers,  particularly  cumulonimbus  clouds, 
appear  at  the  opposite  end  of  the  grey-shade  spectrum  as 
values  in  the  high  50's  or  60's.  Due  to  the  high  albedo  of 
such  clouds,  they  appear  nearly  white  on  the  DMSP  VIS  imagery 

For  DMSP  IR  data,  low  DCV's  and  hence  dark  shades  of  grey 
are  indicative  of  low  temperatures  while  high  DCV's  indicate 
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relatively  high  temperatures.  Because  of  this  convention, 
the  IR  field  is  inverted  prior  to  developing  a positive 
print  (IR  imagery)  so  that  the  cloudy  regions  appear  white, 
and  thus  relate  well  to  that  which  the  human  eye  is 
accustomed  to  viewing  from  the  earth's  surface.  For  both 
VIS  and  IR  DMSP  data  a digital  value  of  63  represents  miss- 
ing data. 

The  data  stored  on  magnetic  tape  were  processed  through 
the  NPS  computer  facility  and  transferred  to  a print-out 
format  for  ease  of  accessing  and  extracting  the  required 
data.  In  order  to  access  the  data,  ship  latitude/longitude 
positions  were  converted  to  I-J  grid  coordinates.  Once  the 
I-J  grid  positions  were  known,  data  were  easily  extracted 
from  the  computer  print-out.  The  formulae  involved  in  the 
conversion  from  latitude/longitude  to  I-J  coordinates  were 
programmed  on  a hand-held  programmable  calculator  and  are 
included  in  Appendix  A.  I-J  grid  positions  were  calculated 
to  the  nearest  hundredth  so  that  an  overlay  subdivided 
accordingly  could  be  used  with  each  I-J  grid  square  on  the 
computer  print-out. 

An  item  of  imp or  tance  involved  the  accuracy  of  extract- 
ing a single  point  value  from  the  available  data.  It  was 
felt  that  an  averaging  scheme  which  took  into  account  sur- 
rounding DCV ' s would  be  more  representative  than  a single 
point  value.  Accordingly,  two  averaging  schemes  were  used 
in  determining  average  DCV 1 s . The  procedures  involved  in 
these  averaging  schemes  are  included  in  Appendix  B.  All 


subsequent  work  with  the  DMSP  VIS  and  IR  data  was  accom- 
plished using  the  "four  by  four"  averaged  parameters. 

The  most  drastically  limiting  factor  involved  in  extract- 
ing the  DMSP  data  proved  to  be  the  time  of  observation. 

Since  the  DMSP  satellite  is  polar  orbiting,  composite  images 
and  digital  data  covering  the  eastern  North  Pacific  Ocean 
area  were  available  only  at  specific  times.  Synoptic  ship 
report  times  of  0000,  0600,  1200  and  1800  GMT  in  many  in- 

stances did  not  coincide  with  times  of  satellite  observation. 
It  was  decided  to  accept  the  DMSP  observation  as  representa- 
tive (in  time)  if  it  was  within  two  hours  of  the  actual 
synoptic  time;  data  outside  the  two-hour  limit  was  rejected. 
With  this  two-hour  criterion,  only  90  VIS  and  196  IR  reports 
could  be  matched  to  the  original  1335  ship  reports  available 
during  the  five-day  period.  These  totals  further  dwindled 
to  69  VIS  and  116  IR  reports  with  the  limitations  imposed 
by  SMS-2  coverage  to  or'iy  the  0000  and  1800  GMT  reports. 

E.  SMS-2  DATA 

SMS-2  digital  data  were  furnished  on  magnetic  tape  by 
the  Space  Science  and  Engineering  Center  (SSEC)  of  the 
University  of  Wisconsin  at  Madison,  Wisconsin. 

The  digital  count  values  range  from  0 to  255  for  both 
VIS  and  IR  data.  In  the  case  of  VIS  DCV's,  low  (high)  num- 
bers correspond  to  surfaces  having  low  (high)  albedo,  while 
for  IR  data  relatively  high  (low)  temperatures  appear  as 
low  (high)  DCV's.  It  would  be  expected,  therefore,  that  a 
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cloud  mass  with  considerable  vertical  development,  such  as 
a cumulonimbus  tower,  would  exhibit  both  VIS  and  IR  values 
at  the  upper  end  of  the  count  value  range.  The  overall 
range  of  DCV's  encountered  in  this  study  was  36-196  for  VIS 
data  and  73-197  for  IR  data. 

The  nominal  accuracy  of  SMS-2  VIS  data  is  one-half  nmi 
by  one-half  nmi  at  the  subsatellite  point  (vicinity  of  Equa- 
tor, 135W) . Actual  values  of  line  and  element  spacing  at 
full  resolution  as  a function  of  latitude/longitude  position 
are  included  in  Fig.  5.  The  variation  of  line  and  element 
values  as  one  proceeds  poleward  and  either  eastward  or  west- 
ward of  the  subsatellite  point  should  be  noted.  The  stated 
SMS -2  accuracy  is  dependent  upon  many  parameters  such  as 
satellite  spin  rate,  satellite  attitude,  orbit  eccentricity, 
camera  geometry,  orbit  precession,  and  height  above  the 
earth's  surface.  Data  resolution  of  two  nmi  by  two  nmi  (i.e. 
1/4  nominal  VIS  resolution)  was  considered  adequate  and 
therefore  used  in  the  processing  of  all  VIS  data,  thus 
greatly  reducing  the  quantity  and  cost  of  data  to  be  pro- 
cessed . 

Reference  has  been  made  to  line  and  element  values  when 
addressing  the  SMS-2  data.  Lines  and  elements  merely  refer 
to  the  manner  in  which  SMS-2  data  are  recorded  on  magnetic 
tape  or  displayed  on  a data  print-out.  Each  VIS  or  IR  DCV 
is  identified  or  indexed  with  an  associated  line  and  element 
number  pair  that  uniquely  locates  the  DCV  on  a projection  of 
the  earth's  surface.  Line  and  element  orientation  are 
depicted  on  a map  of  the  study  area  in  Fig.  6.  It  must  be 
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emphasized,  however,  that  line  and  element  orientation  are 


functions  of  satellite  geometry  and  hardware  and  are,  there- 
by, not  a direct  function  of  parallels  of  latitude  or  merid- 
ians of  longitude  on  the  earth's  surface.  The  two  nmi  by 
two  nmi  data  available  for  this  study  represent  averages  of 
four  data  elements  along  every  fourth  line. 

SMS-2  IR  data  resolution  is  two  nmi  by  four  nmi.  In 
the  format  mode  available  to  the  study,  data  are  spaced  at 
two  nmi  along  each  line  with  each  line  duplicated  once  to 
standardize  the  VIS  and  IR  data  depictions  to  a two  nmi  by 
two  nmi  interval.  Sample  VIS  and  IR  data  are  included  in 
Figs.  7 and  8,  respectively. 

In  addition  to  the  DCV  data  on  magnetic  tape,  the  SSEC 
furnished  a print-out  of  corresponding  line  and  element 
values  at  ten  degree  la ti tude/ long i tude  increments.  Thus, 
an  interpolation  scheme  was  necessary  to  match  the  ship's 
latitude/longitude  positions  (to  the  nearest  tenth)  to  the 
corresponding  line/element  values.  The  solution  appeared 
to  be  in  the  use  of  the  natural  bicubic  spline  curvilinear 
interpolation  scheme  available  at  the  NPS  Computer  Center 
library.  Using  this  routine,  line  and  element  values  could 
be  interpolated  to  the  nearest  whole  degree  of  latitude/ 
longitude  and  linearly  interpolated  thereafter  for  finer 
accuracy . 

Interpolation  using  the  bicubic  spline  was  successfully 
attempted  but  the  processing  time  involved  appeared 
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prohibitive,  especially  since  the  temporal  change  of  all 
satellite  variables  also  necessitated  a new  calculation  of 
line  and  element  grids  for  each  new  synoptic  time.  With 
this  realization,  it  became  evident  that  data  processing 
via  these  methods  was  far  too  time  consuming  and  an  alter- 
nate method  for  handling  the  data  was  sought. 

The  processing  problem  was  resolved  by  extracting  the 
data  required  for  the  study  at  the  SSEC  using  the  Man- 
Computer  Interactive  Data  Access  System  (MCIDAS)  (see 
Appendix  C).  Accordingly,  for  four  days  during  the  latter 
part  of  December  1976,  the  author  processed  data  at  the 
SSEC  using  the  MCIDAS.  The  accessing  of  required  data  via 
line  and  element  number  indices  proved  to  be  a simple  task 
with  the  MCIDAS. 

Having  read  in  the  magnetic  tape,  an  analog  image  of 
the  digital  data  was  generated  on  one  of  the  MCIDAS  ter- 
minals. Once  the  "satellite  picture"  was  generated,  certain 
geographic  features  (i.e.  landmarks)  on  the  earth  were 
identified  and  their  known  latitude/longitude  positions 
were  compared  to  lati tude/ longi tude  positions  that  had  been 
transformed  from  the  satellite  data.  For  the  eastern  North 
Pacific  Ocean,  few  landmarks  were  available  and  were,  of 
necessity,  located  at  the  periphery  of  the  study  area  or 
even  further  inland  on  the  continent.  Landmarks  that  were 
consistently  visible  during  the  5-9  August  time  period  in- 
cluded Isla  Angel  de  la  Guarda  in  the  Gulf  of  California, 
the  Great  Salt  Lake  in  Utah,  Santa  Rosa  Island  off  the 
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southern  California  coast,  and  a noticeable  feature  in  the 
Great  Lakes  region.  During  the  study  period,  none  of  the 
islands  in  the  Aleutian  Chain  could  be  used  as  landmarks 
because  of  the  cloud  cover  consistently  over  that  region. 

For  purposes  of  this  study,  maximum  navigational  errors 
of  10  nmi  were  deemed  acceptable.  When  navigation  informa- 
tion was  checked  for  each  of  the  four  synoptic  times  in  the 
five-day  period,  the  SMS-2  data  at  0600  and  1200  GMT  were 
consistently  outside  acceptable  tolerances  while  0000  and 
1800  GMT  could  be  adjusted  to  within  the  10  nmi  criterion. 

In  addition,  at  synoptic  times  0600  and  1200  GMT  (night-time 
in  the  study  area)  SMS-2  IR  data  only  were  available.  Thus, 
the  deicison  was  made  to  set  aside  the  0600  and  1200  GMT 
data  and  continue  processing  only  the  0000  and  1800  GMT  VIS 
and  IR  data.  At  this  point,  the  conventional  data  set  was 
reduced  to  325  available  ship  reports.  Navigational  adjust- 
ment procedures  as  alluded  to  above  will  not  be  discussed 
herein  due  to  the  detail  involved. 

Once  the  navigation  problem  was  solved,  extraction  of 
necessary  digital  count  value  data  was  a simple  matter. 

Rather  than  extracting  a single  point  value  of  VIS  or  IR 
data,  11  x 11  matrices  of  DCV ' s were  printed  out  for  each 
transient  ship  latitude/longitude  position,  one  matrix  each 
for  the  VIS  and  IR  data,  respectively.  Simple  averages  of 
the  eight  count  values  surrounding  the  center  value  plus 
the  center  point  value  were  also  performed.  These  "nine- 
point"  averages  were  used  in  subsequent  statistical  analyses. 
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Weighted  averages  as  used  with  the  DMSP  data  were  not  com- 
puted due  to  time  limitations  (see  Figs.  7 and  8). 

Another  important  parameter  that  was  derived  based  on 
SMS-2  DCV 1 s was  one  involving  cloud  heights.  This  parameter 
was  available  through  MCIDAS  processing  and  was  extracted 
for  each  ship  latitude/longitude  position.  A brief  discus- 
sion of  the  MCIDAS  Cloud  Height  ( CLDHGT ) program  (Mosher, 
1974),  as  it  applies  to  this  study  is  included  in  Appendix  C 


I 

i 

i 


■: 

I 


i 


I V . DATA  ANALYSIS  AND  RESULTS 


A.  APPROACH 

The  organization  of  data  was  based  on  ground-truth  ship 
report  positions  at  0000  and  1800  GMT  for  the  five-day 
period  5-9  August  1976.  A sample  data  tabulation  is  included 
in  Fig.  2.  . 

The  analysis  of  the  fog/no  fog  data  was  related  to  a 
parameter  developed  in  marine  fog  research  at  the  NPS  (Renard, 
et  al,  1975;  Willms,  1975),  namely  the  fog  duration  category 
(FDC).  Fog  duration  categories  (FDC's)  were  manually  deter- 
mined for  each  of  the  synoptic  ship  reports  based  on  unique 
combinations  of  visibility  plus  present  and  past  weather. 
Willms'  classification  scheme,  based  on  the  visibility  wea- 
ther group  elements,  is  included  in  Table  II.  24  of  Willms' 

39  possible  FDC's  are  represented  in  this  study.  Since  the 
data  set  was  not  significantly  large,  it  w?-  decided  to 
combine  the  FDC's  into  five  groups  for  further  processing 
(see  Table  III).  At  this  point,  the  general  approach  was  to 
examine  the  available  data  with  the  FDC  parameter  acting  as 
the  common  basis  of  comparison.  When  the  SMS-2  data  process- 
ing was  completed,  718  data  points  could  be  matched  to  the 
825  synoptic  ship  reports  available. 

B.  DMSP  DATA  ANALYSIS  AND  RESULTS 

As  previously  indicated,  the  size  of  the  DMSP  data  set 
was  very  small.  Of  718  ground  truth  ship  reports  available, 
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69  VIS  and  116  IR  DMSP  digital  count  values  were  available 
for  analysis.  Using  a histogram  representation,  results 
are  presented  in  Figs.  9 through  12.  For  histograms  that 
were  generated  through  the  NPS  computer  facility,  such  as 
Figs.  9 and  11,  a plotted  empirical  density  function  (F) 
is  included.  A line  through  the  "F"  points  essentially 
gives  a smoothed  profile  of  the  frequency  data.  In  addition, 
histogram  percentage  frequencies  are  indicated  on  the  ordi- 
nate and  the  mean  value  of  each  distribution  is  indicated 
by  a vertical  column  of  "M's".  Table  IV  summarizes  the 
results  concerning  the  distribution  of  DMSP  digital  values. 

The  IR  distribution  (Fig.  11)  appears  to  approach  a 
normal  distribution  to  a much  greater  degree  than  does  the 
VIS  counterpart  (Fig.  9),  but  lack  of  sufficient  data  makes 
such  a conclusion  very  tenuous.  The  range  of  VIS  DCV ' s 
(26  - 60)  is  somewhat  greater  than  that  for  IR  DCV ' s (19  - 
46)  with  the  distribution  more  erratic  for  the  former. 

Similar  observations  are  appropriate  when  considering  each 
of  the  FDC  1 s in  Figs.  10  and  12  and  Table  IV. 

C.  SMS-2  DATA  ANALYSIS  AND  RESULTS 

1 . SMS-2  VIS  and  IR  DCV  Analysis 

- Histograms  and  Scatter  Diagrams 

The  SMS-2  data  set  proved  to  be  more  informational 
than  the  DMSP  data,  due  in  part  to  the  significant  increase 
in  size  of  the  data  set.  It  must  be  emphasized,  however, 
that  when  dealing  with  a sample  size  numbering  718,  sta- 
tistical findings  are  necessarily  more  indicative  or 


suggestive  than  truly  definitive  as  compared  with  a sample 
size  numbering  say,  10,000*  With  this  caution  in  mind, 
results  of  the  analysis  of  the  SMS-2  data  follow. 

As  in  the  case  of  the  DMSP  data,  all  SMS-2  data, 

VIS  and  IR  separately,  were  organized  into  the  same  five 
groups  of  FDC's.  Figs.  13  and  19  depict  the  distributions 
of  VIS  and  IR  average  DCV's  respectively,  for  all  data. 
Average  count  values  were  selected  as  more  representative 
than  point  values.  However,  average  count  value  distribu- 
tions were  not  significantly  different  from  their  point- 
value  distribution  counterparts.  As  had  been  observed  with 
the  DMSP  data,  a difference  in  range  interval  was  noted  in 
VIS  as  opposed  to  IR  data.  VIS  data  encompassed  159 
possible  values  while  IR  data  spanned  123  possible  digital 
count  values.  In  fact,  for  all  groups,  except  FDC  GR  2, 
the  IR  DC V ranges  are  less  than  those  observed  for  their 
VIS  counterparts.  It  was  also  interesting  to  note  that  the 
absolute  minimum/maximum  count  values  encountered  were  36- 
198  for  VIS  and  73-196  for  IR;  possible  SMS-2  digital  count 
values  are  advertised  to  range  from  0-255. 

After  the  entire  data  set  had  been  scrutinized  on 
the  initial  two  plots,  additional  histograms  segregating 
the  data  into  the  FDC  groups  were  constructed  and  appear 
as  Figs.  14  through  18  and  20  through  24.  Additionally, 
Table  V summarizes  the  essential  statistics  from  all  of  the 
SMS-2  histogram  data.  The  following  additional  observations 
were  made  based  on  the  histograms: 
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(1)  Fog  reports  seldom  occur  in  the  VIS  count 
value  range  of  approximately  36-72  (i.e.  low 
albedo)  (see  Figs.  14  through  17). 

(2)  Similarly,  fog  reports  seldom  occur  in  the 

IR  count  value  range  of  73-89  (i.e.  relative- 

ly high  temperatures )( see  Figs.  20  through  23). 

(3)  Limiting  values,  similar  to  (1)  and  (2)  above, 
at  the  upper  end  of  the  VIS  (high  albedo)  and 
IR  (low  temperature)  count  value  range  are 
not  evident. 

(4)  Mean  DCV 1 s differ  little  in  all  categories 
associated  with  fog  for  both  VIS  and  IR  data 
(Table  V ) . 

The  presentation  of  any  one  parameter  as  a function 
of  FDC  groups  was,  in  essence,  a one-dimensional  view  of 
the  data.  To  further  pursue  a solution  to  the  basic  objec- 
tive of  the  research,  a multi-dimensional  data  depiction 
method  was  desirable.  Accordingly,  two-dimensional  plots 
were  constructed  with  VIS  and  IR  DCV's  as  the  two  axes. 

For  each  ground-truth  ship  latitude/longitude  position,  the 
VIS  average  DCV  was  plotted  versus  the  IR  average  DCV.  Five 
plots  were  generated  (Figs.  25  through  29)  which  segregated 
the  data  into  the  five  FDC  groupings  previously  defined. 
Information  depicted  on  these  plots  was  identical  to  the 
histogram  data  but  allowed  VIS  and  IR  data  to  be  "spatially" 
correlated. 

Objective  analysis  techniques  employed  up  to  this 
point  had  still  not  established  any  limited  range  of  either 


VIS  or  IR  count  values  for  fog  occurrence.  In  an  effort  to 
establish  an  upper  cut-off  range  of  VIS/IR  values  for  fog 
occurrence  (to  supplement  the  limiting  low  values  already 
determined)  a subjective  analysis  technique  was  selectively 
employed  using  both  VIS  and  IR  satellite  imagery.  The 
principle  involved  examination  of  all  data  locations,  except 
those  associated  with  FDC  180  (no  fog)  on  the  VIS/IR  imagery 
in  an  effort  to  segregate  out  those  reports  in  which  middle 
and/or  high  clouds  were  in  evidence.  This  step  did  not 
negate  surface-based  observations  of  fog  but  merely  identi- 
fied those  ship  locations  where  the  satellite  was  seeing  a 
low  cloud  type  (as  stratus)  only. 

The  accomplishment  of  this  subjective  "sanitization" 
proved  the  necessary  breakthrough  in  establishing  reasonable 
cutoff  boundaries  for  the  IR  count  values.  Figs.  30  through 
33  depict  the  "sanitized"  count-value  distributions  for 
the  fog-associated  FDC  groups.  A comparison  of  these  fig- 
ures with  the  immediately  preceding  set  (Figs.  25  through 
28)  graphically  illustrates  the  significant  decrease  in  IR 
count  value  range.  Table  VI  summarizes  SMS-2  count  value 
parameters  before/after  the  subjective  sanitization  pro- 
cedure. It  is  important  to  note  the  relatively  small  change 
in  the  VIS  range  of  count  values  as  compared  to  the  IR 
range.  Results  observed  are  reinforced  to  some  degree  by 
the  findings  of  Hale  and  Renard  in  their  1975  study  which 
utilized  NOAA-2  digitized  data  in  diagnosing  marine  fog. 

A large  VIS  range  of  possible  DCV's  and  a much  more  limited 
range  of  IR  count  values  were  reported  in  their  study. 
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During  the  selective  sanitization  of  the  ship  lati- 
tude/longitude positions,  it  was  observed  that  a certain 
number  of  reports  when  viewed  on  the  VIS/IR  imagery  appeared 
to  be  in  areas  of  large  gradients  of  grey  shade  values,  and 
hence  the  interpretation  became  critically  dependent  on 
gridding  and  navigational  accuracies  as  well  as  on  the  obser- 
vation itself.  Two  examples  are  included  here  to  illustrate 
the  high  observed  correlation  between  imagery  and  associated 
digital  data. 

The  first  example,  depicted  in  Fig.  34,  is  from  0000 
GMT  data  on  5 August  1976.  The  VIS  and  IR  imagery  at  the 
ship  latitude/longitude  position  indicated  a uniform  layer 
of  low  clouds  with  several  localized  middle  and/or  high  cloud 
patches  surrounding  the  ship  position.  The  VIS  data  matrix 
reflects  this  condition  quite  accurately  in  the  strong  count 
value  gradients  extant  across  the  data  matrix.  The  IR  data 
matrix,  due  primarily  to  resolution  differences, does  not 
exhibit  the  gradient  as  well  except  at  the  right  edge  of  the 
matrix.  To  further  confirm  the  satellite  data,  the  synoptic 
ship  report  indicated  f ractostratus  reported  for  low  cloud 
type  and  double-layered  altocumulus  and  dense  cirrus  in 
patches  for  the  middle  and  high  cloud  types  respectively. 

The  dimensions  of  the  SMS-2  VIS/IR  data  matrices  (Fig.  34) 
are  approximately  17  nmi  (approximately  east/west)  by  32 
miles  (approximately  north/south). 

The  VIS  and  IR  data  depicted  in  Fig.  35  correspond 
to  an  area  on  the  VIS  imagery  that  constituted  a "clear  slot" 
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crossing  precisely  at  the  ship  latitude/longitude  position 
indicated.  The  strong  count  value  gradients  observed  on 
either  side  of  the  slot  indicate  the  presence  of  highly  re- 
flective clouds.  The  IR  data  does  not  reflect  this  "clear 
slot”  probably  due  to  the  previously  mentioned  resolution 
differences  and  the  resultant  "smearing"  of  the  digital 
data.  VIS/IR  data  matrix  dimensions  at  the  geographical 
position  are  approximately  20  nmi  (approximately  east/west) 
by  32  nmi  (approximately  north/south). 

2 . SMS-2  Cloud  Height  Parameters 

Since  the  histogram  approach  of  data  presentation 
had  proved  its  utility,  it  was  employed  again  in  an  effort 
to  glean  some  meaning  from  the  SMS-2  cloud  height  parameters 
which  had  been  generated  using  a program  developed  by 
Mosher  at  the  University  of  Wisconsin  (Mosher,  1974).  As 
discussed  in  Appendix  C,  the  parameter  of  interest  was  the 
derived  height  of  the  cloud  base.  Included  in  Figs.  36 
through  41  are  histograms  of  cloud  base  height  in  meters 
(abscissa)  as  a function  of  the  various  FDC  groups  previously 
defined.  The  CLDHGT  program  had  not  been  developed  for 
application  in  fog  research  and  hence  was  not  specifically 
"tuned"  to  necessarily  provide  cloud  base  heights  near  zero 
as  one  observes  with  fog.  Some  of  the  assumptions  inherent 
in  the  CLDHGT  program  such  as  the  development  and  tailoring 
of  the  model  to  cumulus-type  clouds  obviously  would  affect 
generated  results.  It  was  not  a great  surprise  then,  that 
negative  heights  were  generated  using  the  CLDHGT  program. 
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Perhaps  the  most  meaningful  inference  drawn  from 
the  cloud  height  histograms  is  the  near  Gaussian  distribu- 
tion (without  regard  to  negative  cloud  heights)  observed 
particularly  in  Figs.  36,  37  and  39.  It  should  be  noted 

that  mean  cloud  base  heights  calculated  for  each  of  the 
fog  groups  were,  in  each  instance,  less  than  the  mean 
height  calculated  for  the  FDC  180  (no  fog)  group.  Table 
VII  summarizes  essential  results  from  the  CLDHGT  program. 

The  full  impact  that  a satellite-generated  cloud 
height  program  might  have  in  achieving  the  goals  stated 
for  this  study  cannot  be  answered  by  this  research.  Further 
investigation,  refinement,  and  possible  "tuning"  of  the 
CLDHGT  program  for  this  application  to  marine  fog  should 
be  accomplished. 

3 . Other  Statistical  Relationships 

To  supplement  the  fog  study  investigation  without 
specific  application  of  the  digital  satellite  data,  weather 
parameters  available  in  the  standard  synoptic  ship  report 
were  also  statistically  examined.  Using  the  histogram 
approach  again,  air  temperature,  sea  surface  temperature, 
and  surface  wind  direction  were  depicted  as  functions  of 
the  various  FDC  groups.  The  parameter  that  proved  most  in- 
teresting was  the  distribution  of  the  surface  wind  as  a 
function  of  FDC.  The  surface  wind  distribution  for  those 
data  in  category  FDC  GR  1 will  be  specifically  addressed 
here . 

In  examining  the  histogram  initially  (Fig.  42), 
some  two-thirds  of  the  wind  reports  fell  between  directions 
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150  and  270  deg.  To  depict  the  surface  wind  data  in  a 
second  fashion  and  to  also  indirectly  involve  the  VIS  and 
IR  DCV ' s , wind  directions  were  attached  to  each  of  the 
data  points  depicted  in  Fig.  31.  To  establish  mean  wind 
directions,  averages  by  column  and  row  as  portrayed  in 
Fig.  43  were  performed.  For  the  averages  given  in  which 
five  or  more  wind  reports  were  considered  in  the  average, 
the  mean  wind  directions  were  as  follows;  160  through 
262  along  the  VIS  axis;  171  through  226  along  the  IR  axis. 
Thus,  the  surface  winds  with  southerly  components  are  typi- 
cally related  to  heavy  fog  conditions  over  the  eastern 
North  Pacific  Ocean  during  the  August  period  studied. 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  principal  objective  of  this  study  was  to  ascertain 
the  value  and  applicability  of  DMSP  and  SMS-2  satellite  data 
in  the  diagnosis  of  marine  fog  over  the  open  ocean.  Through 
the  interpretation  and  processing  of  digital  data  from  both 
satellite  systems  with  conventional  surface  ship  reports  as 
ground-truth  verification,  optimal  ranges  of  digital  count 
values  related  to  fcg  occurrence  for  both  the  VIS  and  IR 
data  were  determined  i-'ble  Vi).  A secondary  objective  of 
the  study,  the  establishment  of  usable  statistical  correla- 
tions of  satellite  and  conventionally-reported  weather  para- 
meters, was  only  partially  achieved  through  the  relationship 
between  observed  fog  reports,  associated  surface  wind,  and 
SMS -2  VIS  and  IR  DCV's.  Further  statistical  analyses  were 
not  performed  using  these  data  due  to  time  limitations. 

DMSP  data  proved  insufficient  in  quantity  for  truly 
meaningful  input  to  the  study  as  a whole.  Any  polar  orbit- 
ing satellite  data  utilized  in  a study  such  as  this  one 
must  be  available  over  a much  longer  time  period  so  that  a 
representative  number  of  data  points  may  be  made  available 
for  investigation.  Processing  techniques  employed  herein, 
however,  may  be  useful  for  future  studies  involving  DMSP 
data . 

SMS-2  recorded  data  processed  using  the  MCIDAS,  pro- 
vided the  major  data  source  for  the  study.  The  high  quality 
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numerical  products  and  imagery  available  from  the  SMS-2 
system  provided  the  necessary  base  for  meaningful  data  pro- 
cessing. Even  though  the  final  sample  size  for  the  study 
was  relatively  small,  the  processing  accuracies  attained 
throughout  the  investigation  were  noteworthy  and  provided 
very  high  quality  data  for  further  analysis.  The  utility 
of  a derived  cloud  height  parameter  that  has  been  finely 
tuned  to  account  for  various  cloud  types  and  varied  atmos- 
pheric conditions  is  apparent.  Although  the  MCIDAS  CLDHGT 
parameters  derived  for  this  study  did  not  provide  signifi- 
cant insight  into  fog  diagnosis,  it  is  felt  that  refinement 
and  development  of  this  program  will  provide  valuable  input 
toward  future  studies.  Perhaps  the  CLDHGT  program  will 
also  be  able  to  address  still  another  problem,  namely,  the 
segregation  and  discrimination  of  fog  from  stratiform  clouds 
without  associated  fog. 

The  following  recommendations  are  offered  for  future 
studies : 

(1)  Conduct  a similar  study  for  a longer  time  period 
using  both  DMSP  and  SMS-2  data. 

(2)  Work  toward  increased  automation  of  processing 
techniques  herein  outlined  using  existing  computer 
facilities  . 

(3)  Investigate  the  feasibility  of  modifying  existing 
MCIDAS  navigational  and  CLDHGT  programs  for  use  on 
the  NPS  computer. 
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Further  investigate  the  MCIDAS  CLDHGT  program  as 


developed  by  Mosher  (1974)  in  the  fog  diagnosis 
application . 

Further  examine  statistical  relationships  be- 
tween conventional  synoptic  ship  report  data  and 
VIS/IR  count  value  distributions. 

Consider  future  fog  research  using  the  MCIDAS 
facility  at  the  SSEC,  University  of  Wisconsin. 
Investigate  a scheme  which  would  involve  the 
assignment  of  "weighting  factors"  to  selected 
VIS/IR  digital  count  value  parameters  and  when 
used  in  combination  with  climatology,  moisture 
content,  temperature,  wind,  and  other  fog- 
related  parameters,  arrive  at  a credible  fog 
probability  analysis  (Hale  and  Renard,  1975)  . 


A 
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APPENDIX  A 


Conversion  of  latitude/longitude  position  to  I-J  grid  co- 
ordinates (FNWC  Subroutine  Library). 

The  following  expressions  were  programmed  for  use  with 

a programmable  calculator: 


1 = *p  + t<Re>  (i  + sin  (cos ( 350-X) ) ] 
J = Jp  + t(Re)  (l  r-'sTn  ^ (sin(350-X))] 


where:  1^  = Pole  position  = 31 

Jp  = Pole  position  = 31 
<p  = Latitude,  in  degrees 
A = Longitude,  in  degrees 

Re  = Distance  from  pole  to  equator  in  grid 
mesh  lengths  = 31.205. 
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APPENDIX  B 


Averaging  schemes  for  DMSP  digital  data  - 2x2  and  4x4. 

Given : 

(1)  Set  of  I-J  grids  computer  produced,  each  consisting 

of  64  by  64  (4096  total)  digital  count  values. 

(2)  Latitude/longitude  positions  converted  to  I-J  grid 
coordinates  (see  Appendix  A). 

(3)  Overlay,  subdivided  to  nearest  hundredth,  to  extract 
values  from  printed  grid. 

Procedure : 

(1)  Enter  I-J  grid  (with  overlay)  and  select  the  two- 
digit  value  closest  to  point  on  the  overlay. 

(2)  Extract  surrounding  four  points  (including  point 
value  just  selected) ; sum  points  and  divide  by 
four  for  2x2  average. 

(3)  For  calculation  of  the  4x4  average  a weighting 
scheme  is  used, based  on  radial  distance  from  the 
center  of  the  grid.  I-J  coordinates  entered  a-e 
assumed  to  be  at  grid  center  to  accomplish  this 
average . 
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Points  a,  b,  c refer  to  digital  count  values  as  located  on 
64x64  grid. 


t 


Assume:  Distance  between  elements  is  equal  to  two  units 

for  purposes  of  deriving  weighting  factors  (W  ) , (i.e. 

distance  a to  b =2  units). 
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Sample  Averaging  Procedure: 


= Exact 
point  value 

2x2  average 

4x4  average 
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location  of  I-J  coordinate  (to  nearest  hundredth) 
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APPENDIX  C 


I 


The  Man-Computer  Interactive  Data  Access  System  (MCIDAS) 

The  Man-Computer  Interactive  Data  Access  System  (MCIDAS) 
located  at  the  Space  Science  and  Engineering  Center  at  the 
University  of  Wisconsin,  Madison,  Wisconsin,  was  developed 
for  processing  two  dimensional  image  data  from  high  volume 
data  sources  such  as  the  Synchronous  Meteorological  Satellite 
(SMS-2)  which  transmits  approximately  75  billion  bits  of 
image  data  each  day  (Smith,  1975).  The  characteristics  that 
are  required  for  a system  to  process  SMS-2  data  and  which 
were  therefore  built  into  the  MCIDAS  include: 

"a.  Best  possible  navigation 

b.  Preservation  of  inherent  resolution  and  geometry 

c.  Speed  and  precision  of  automated  techniques 

d.  Human  interaction  only  for  decision-making  but  not 
for  data  handling  or  measurement 

e.  Combination  of  digital  and  analog  hardware 

f.  Real  or  near-real  time  processing." 

(Chatters  and  Suomi,  1975) 

Some  of  the  capabilities  of  the  MCIDAS  include  the 
following : 

(1)  Digital  data  may  be  transformed  into  analog  picture 
format  display  at  MCIDAS  terminals. 

(2)  Images  may  be  time  sequenced  on  a terminal;  a use- 
ful capability  in  tracking  clouds. 

(3)  Inherent  in  the  system's  operation  is  an  interactive 
capability  in  which  the  human  operator  always 
occupies  a decision-making  position. 
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(4)  Images  (pictures)  on  a MCIDAS  terminal  may  be  com- 
bined or  overlaid  as  well  as  color  enhanced. 

Unique  features  which  are  intrinsic  to  the  MCIDAS  include 
both  a unique  navigational  model  and  the  WINDCO/CLDHGT  system. 
The  navigational  model,  although  extremely  complex,  allows 
rapid  transformation  from  satellite  coordinates  to  true 
earth  coordinates.  The  WINDCO/CLDHGT  system  was  developed 
for  use  in  cloud-tracking . 

For  this  study,  the  CLDHGT  sub-system  provided  a unique 
parameter:  cloud  base  height  derived  from  satellite  digital 

data.  The  basic  steps  involved  in  the  program  include  the 
following : 

(1)  Using  visible  brightness  values  the  optical  thick- 
ness of  a cloud  is  determined. 

(2)  From  the  optical  thickness,  the  IR  emissivity  of 
the  cloud  is  computed. 

(3)  The  true  temperature  of  the  cloud  is  then  computed. 

(4)  From  this  temperature,  and  climatological  soundings 
indexed  by  season  and  geographical  position,  the 
height  of  the  cloud  top  is  determined. 

(5)  Cloud  thickness  can  be  computed  from  the  optical 
thickness  and  an  extinction  coefficient  (Mosher, 

1974)  . 

Cloud  thickness  was  subtracted  from  cloud  top  height  to  pro- 
duce a cloud  base  height  for  each  ground  truth  ship  latitude/ 
longitude  position.  Since  the  program  was  originally  de- 
veloped and  modeled  for  cumulus-type  clouds,  a tentative  cor- 
rection factor  was  introduced  by  Mr.  Mosher  for  the  purpose 
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Figure  2.  Sample  tabulated  data,  0000  GMT,  0 August  1976 
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Figure  7.  Sample  SMS-2  grid  of  VIS  digital  count 

values  (DCV's);  center  value  (circled)  is 
at  ship  latitude/longitude  position; 
boxed  area  includes  nine  points  used 
to  obtain  average  DCV . 
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Figure  3, 


Sample  SMS-2  grid  of  IR  digital  count 
values  (DCV's) j center  value  (circled)  is 
at  ship  latitude/longitude  position; 
boxed  area  includes  nine  points  used 
to  obtain  average  DCV. 
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Figure  9.  Frequency  distribution  of  DMSP  VIS  average  digital  count  values  (DCV's) 
for  all  fog  duration  categories  (see  Tables  II  and  III);  30-60N, 
eastern  North  Pacific  Ocean,  b-9  August  1976;  (1)  frequencies  and 

sample  size  along  abscissa  (top),  (2)  percentage  frequencies  along 
ordinate,  (3)  M = mean  DCV , ( 4 ) F - empirical  density  function. 


Figure  Jl.  Frequency  distribution  of  DMSP  Ifi  average  digital  count  values  (DCV's) 
for  all  fog  duration  categories  (see  Tables  II  and  III);  30-60N, 
eastern  North  Pacific  Ocean,  5-9  August  1976;  (1)  frequencies  and 

sample  size  along  abscissa  (top),  (2)  percentage  frequencies  along 
ordinate,  (3)  M = mean  DCV , (4)  V - empirical  density  function. 
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Figure  13.  frequency  distribution  of  SMS-2  VIS  average  digital  count  values  (DCV's) 
for  all  fog  duration  categories  (see  Tables  II  and  III);  30-60N,  eastern 
North  Pacific  Ocean,  5-9  August  1976;  (1)  frequencies  and  sample  size 

along  abscissa  (top),  (2)  percentage  frequencies  along  ordinate,  (3)  M = 
mean  DCV,  (4)  P = empirical  density  function. 
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Figure  16.  Frequency  distribution  of  SMS-2  VIS  average  digital  count  values  (DCV's) 
for  fog  duration  category  group  3 (see  Tables  II  and  III);  30-60N, 
eastern  North  Pacific  Ocean,  5-9  August  1976;  (1)  frequencies  and  sample 

size  along  abscissa  (top),  (2)  percentage  frequencies  along  ordinate, 

(3)  M = mean  DCV,  (4)  F = empirical  density  function. 
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Frequency  distribution  of  SMS-2  VIS  average  digital  count  values  (DCV's) 
for  fog  duration  category  180  (see  Tables  II  and  III);  30-60N,  eastern 
North  Pacific  Ocean,  5-9  August  1976;  (1)  frequencies  and  sample  size 

along  abscissa  (top),  (2)  percentage  frequencies  along  ordinate,  (3)  M = 
mean  DCV,  (4)  F = empirical  density  function. 
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19.  Frequency  distribution  of  SMS-2  Il<  average  digital  count  values  (DCV's) 
for  all  fog  duration  categories  (see  Tables  II  and  III);  30-60N, 
eastern  North  Pacific  Ocean,  5-9  August  1976;  (1)  frequencies  and 

sample  size  along  abscissa  (top),  (2)  percentage  frequencies  along 
ordinate,  (3)  M = mean  DCV,  (4)  F = empirical  density  function. 


figure  20.  Frequency  distribution  of  SMS-2  Il(  average  digital  count  values  (DCV's) 
for  fog  duration  category  group  1 (see  Tables  II  and  III);  30-60N, 
eastern  North  Pacific  Ocean,  5-9  August  1976;  (1)  frequencies  and 

sample  size  along  abscissa  (top),  (2)  percentage  frequencies  along 
ordinate,  (3)  M = mean  DCV,  (4)  F = empirical  density  function. 


Figure  22.  Frequency  distribution  of  SMS-2  IR  average  digital  count  values  (DCV's) 
for  fog  duration  category  group  3 (see  Tables  II  and  III);  30-60N, 
eastern  North  Pacific  Ocean,  5-9  August  1976;  (1)  frequencies  and 

sample  size  along  abscissa  (top),  (2)  percentage  frequencies  along 
ordinate,  (3)  M = mean  DCV,  (4)  F - empirical  density  function. 


figure  23.  Frequency  distribution  of  SMS-2  IR  average  digital  count  values  (DCV's) 
for  fog  duration  category  175  (see  Tables  II  and  III);  30-60N,  eastern 
North  Pacific  Ocean,  5-9  August  1976;  (1)  frequencies  and  sample  size 

along  abscissa  (top),  (2)  percentage  frequencies  along  ordinate, 

(3)  M = mean  DOV,  (4)  F = empirical  density  function. 


Figure  24.  Frequency  distribution  of  SMS-2  IK  average  digital  count  values  (DCV's) 
for  fog  duration  category  180  (see  Tables  II  and  III);  30-60N,  eastern 
Nortli  Pacific  Ocean,  5-9  August  1976;  (1)  frequencies  and  sample  size 

along  abscissa  (top),  (2)  percentage  frequencies  along  ordinate, 

(3)  M = mean  L)CV , (4)  F = empirical  density  function. 
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Figure  31.  Scatter  diagram  of  SMS-2  VIS  vs.  IR  average  d 
for  fog  duration  category  group  2 (see  Tables 
high  cloud  cases  removed;  30-60N,  eastern  Nor 


Scatter  diagram  of  SMS-2  VIS  vs.  IR  average  digital  count  val 
for  fog  duration  category  175  (see  Tables  II  and  III);  middle 
high  cloud  cases  removed;  30-60N,  eastern  North  Pacific  Ocean 
August  1976. 
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Figure  34.  VIS  (top)  and  IR  bottom)  digital  count  value 
matrices  for  surface  ship  location  43. 3N 
143. 1W  (center  of  data  matrix),  0000  GMT, 

5 August  1976. 
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Figure  36.  Frequency  distribution  of  SMS-2  cloud  base  height  (in  meters)  foi 
fog  duration  categories  (see  Tables  II  and  III);  30-60N,  eastern 
North  Pacific  Ocean,  5-9  August  1976;  (1)  frequencies  and  sample 

along  abscissa  (top),  (2)  percentage  frequencies  along  ordinate, 
(3)  M = mean  l)CV,  (4)  F = empirical  density  function. 
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along  abscissa  (top),  (2)  percentage  frequencies  along  ordinate 
(3)  M = mean  DCV,  (4)  F = empirical  density  function. 
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duration  category  175  (see  Tables  II  and  III);  30-60N,  eastern  North 
Pacific  Ocean,  5-9  August  1976;  (1)  frequencies  and  sample  size  along 

abscissa  (top),  (2)  percentage  frequencies  along  ordinate,  (3)  M = 
mean  DCV,  (4)  F - empirical  density  function. 
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Figure  41.  Frequency  distribution  of  SMS-2  cloud  base  height  (in  meters)  for  fog 
duration  category  180  (see  Tables  II  and  III);  30-60N,  eastern  North 
Pacific  Ocean,  5-9  August  1976;  (1)  frequencies  and  sample  size  along 

abscissa  (top),  (2)  percentage  frequencies  along  ordinate,  (3)  M = 
mean  DCV,  (4)  F = empirical  density  function. 
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Figure  42.  Frequency  distribution  of  surface  winds  for  synoptic  ship  reports; 

fog  duration  category  group  1 (see  Tables  II  and  III);  30-60N, 
eastern  North  Pacific  Ocean,  5-9  August  1976;  (1)  frequencies  and 

sample  size  along  abscissa  (top),  (2)  percentage  frequencies  along 
ordinate,  (3)  M = mean  DCV , (4)  F = empirical  density  function. 
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TABLE  I 

DMSP  data  interval  as  a function  of  latitude 


Latitude 


I-J  grid  size  (nmi) 


Data  interval  (nmi ) 


Equa  tor 

110.4 

by 

110.4 

1.73 

3 ON 

165.6 

by 

165.6 

2 .59 

45N 

188.5 

by 

188 . 5 

2.95 

6 ON 

206 

by 

206 

3.22 

r 


88 


TABLE  II.  Marine-fog  duration  categories  as  dependent  on 
the  visibility-weather  group  elements  of  the  primary  synoptic 
report.  Hours  of  fog  assigned  are  derived  from  an  analysis 
of  summer-season  North  Pacific  Ocean  Weather  Station  (P,Q,S) 
weather  log  data  (after  Willms,  1975). 


Fog 

Duration 

Categories 


Coded 

Visibility 


Coded  Coded 

Present  Past 

Weather  Weather 


* 

* * 
H 


10 

★ 

40-49 

4 

5.5 

15 

* 

40-49 

5 

4.3 

20 

* 

10 

4 

5.1 

22 

* 

11,12 

4 

3.0 

24 

★ 

10 

5 

3.0 

26 

* 

11, 12 

5 

2.0 

30 

90-93 

* * 

4 

4.7 

31 

90-93 

H 

4 

3 . 0 

32 

94 

* * 

4 

3.8 

33 

94 

H 

4 

3 . 0 

35 

90-93 

* * 

5 

3.8 

36 

90-93 

H 

5 

2 . 0 

37 

94 

★ * 

5 

1.2 

38 

94 

H 

5 

1 . 0 

40 

* 

28 

4- 

5 . 3 

45 

★ 

28 

5,6 

3 . 3 

50 

* 

40-49 

6 

3 . 6 

60 

★ 

10 

6 

2 . 6 

65 

★ 

11, 12 

6 

1.0 

70 

90-93 

* * 

6 

2.9 

71 

90-93 

H 

6 

1.0 

72 

94 

★ ★ 

6 

1.0 

73 

94 

H 

6 

1.0 

80 

★ 

40-49 

2 

4.0 

90 

* 

10 

2 

3 . 5 

95 

★ 

11,12 

2 

1 . 0 

100 

90-93 

* * 

2 

2 . 5 

102 

94 

* * 

2 

1.3 

110 

★ 

28 

2 

4.6 

120 

★ 

* * 

4 

3 . 1 

130 

★ 

40-49 

0,1,3, 

7-9 

3.6 

140 

* 

10 

0,1,3, 

7-9 

2.6 

145 

★ 

11,12 

0,1,3, 

7-9 

1.0 

150 

90-93 

* * 

0,1,3, 

7-9 

1.0 

152 

94 

* * 

0,1,3, 

7-9 

1.0 

160 

* 

28 

0,1,3, 

7-9 

3 . 3 

175 

95-96 

* * 

0-3, 

5-9 

0 . 6 

180 

All  reports 

not 

fitting 

above  categories 

Deno  tes 

any  visibility 

allowed,  i 

.e.  90-99 

Deno  tes 

any  present  weather  code 

except  10 

,11,12 

,28 

or 

heavy 

present  weather 

codes 

Denotes 

heavy  present 

weather  codes  only: 

30-39 , 

62 

67, 

69, 

72-75,  81,  82, 

84 

, 86,  88 

, 90-99 

89 


TABLE  III.  Fog  duration  category  (FDC)  groupings. 

Assumed 


Group 

Weather  Condition 

FDC  ' s 

included 

FDC  Group 

1 

heavy  fog  likely 

10,15, 

30,35,50, 

(FDC  GR  1) 

at  observation 
time 

70,80, 

100,130 

FDC  Group 

2 

light  fog  likely 

20,24, 

32,37,60, 

(FDC  GR  2) 

at  observation 
time 

90,140 

FDC  Group 

3 

fog  likely  but  not 

22,26, 

40,73,110 

(FDC  GR  3) 

necessarily  at  ob- 
servation time 

120 

FDC  175 

visibility  qualifies 
for  light  fog,  but 
fog  not  reported 

175 

TABLE  V.  SMS-2  histogram  statistics. 


Range  of  Total  no.  Mean 

Histogram Figure  Count  Values  of  reports  value 


VIS 

All 

FDC 

13 

39 

198 

715 

121.7 

VIS 

FDC 

GR  1 

14 

74 

- 

184 

131 

135.6 

VIS 

FDC 

GR  2 

15 

92 

- 

160 

24 

135.9 

VIS 

FDC 

GR  3 

16 

62 

- 

184 

32 

134.2 

VIS 

FDC 

175 

17 

47 

- 

198 

41 

142.8 

VIS 

FDC 

180 

18 

39 

- 

185 

457 

114.0 

IR 

All 

FDC 

19 

73 

- 

196 

714 

115.6 

IR 

FDC 

GR  1 

20 

89 

- 

195 

131 

121.3 

IR 

FDC 

GR  2 

21 

95 

- 

191 

24 

124.6 

IR 

FDC 

GR  3 

22 

96 

- 

180 

32 

118.8 

IR 

FDC 

175 

23 

76 

- 

190 

41 

129.6 

IR 

FDC 

180 

24 

73 

— 

196 

456 

112.1 

i 

! 


92 


TABLE  VI.  SMS-2  VIS/IR  count  value  range  before/after 
removal  of  middle/high  clouds. 


No . of 

reports  after  Digital  count 

apparent  middle  value  range 


FDC 

Group 

NO  . of 
reports 

and  high  cloud 
cases  removed 

VIS 

before  after 

IR 

before  after 

FDC  GR 

1 

131 

74 

74-184 

78-168 

89-195 

89-130 

FDC  GR 

2 

24 

16 

92-160 

92-160 

95-191 

95-132 

FDC  GR 

3 

32 

23 

62-184 

62-167 

96-180 

96-120 

FDC  175 

41 

18 

47-198 

47-168 

76-190 

76-126 

93 


TABLE  VII 


SMS-2  cloud  base  height  statistics 


Histogram 

Figure 

Mean  height  (m) 

Total  nc 

All 

FDC 

36 

651.5 

7] 

FDC 

GR  1 

37 

462.9 

12 

FDC 

GR  2 

38 

-121.8 

2 

FDC 

GR  3 

39 

287.6 

- 

FDC 

175 

40 

-842.6 

A 

FDC 

180 

41 

915.2 

45 

94 
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